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Abstract. We present an experimental investigation of the agglomeration of microbubbles into a 2D mi- 
crofoam and its flow in a rectangular microchannel. Using a flow-focusing method, we produce the foam 
tn situ on a microfluidic chip for a large range of liquid fractions, down to a few percent in liquid. We 
can monitor the transition from separated bubbles to the desired microfoam, in which bubbles are closely 
packed and separated by thin fllms. We flnd that bubble formation frequency is limited by the liquid flow 
rate, whatever the gas pressure. The formation frequency creates a modulation of the foam flow, rapidly 
damped along the channel. The average foam flow rate depends non-linearly on the applied gas pressure, 
displaying a threshold pressure due to capillarity. Strong discontinuities in the flow rate appear when the 
number of bubbles in the channel width changes, reflecting the discrete nature of the foam topology. We 
also produce an ultra flat foam, reducing the channel height from 250 /im to 8 ^m, resulting in a height 
to diameter ration of 0.02; we notice a marked change in bubble shape during the flow. 
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Two phase microflows like microemulsions, microbubbles 
and microdrops presently attract considerable attention 
PEIEHE]- To these we would like to add microfoams. 
Their application in a lab-on-a-chip context provides the 
possibility for the efficient handling of series of gas pock- 
ets, and allows to create microchemical reactors that are 
both very rapid and highly parallelized. Specifically, the 
gas-liquid interface of microfoams provides a transport lo- 
cation for amphiphilic molecules, with a hydrophilic head 
and a hydrophobic tail. A decrease in size increases the 
surface to volume ratio; hence, microfoams could be used 
as an efficient carrier for proteins or lipids at high concen- 
tration. 

Microfoams offer advantages compared to foams at 
larger scales for the study of foam properties as micro- 
foams are very stable and well controlled. First, the ab- 
sence of vertical drainage on the small length scales of 
a microfluidic system creates liquid profiles in the foam 
films that are constant over time and do not show an 
asymmetry due to gravity. Another advantage is that be- 
cause of the low Reynolds numbers involved, the amount 
of gas produced during bubbling is very stable with nearly 
monodisperse bubble volume distribution. The amount of 
liquid (liquid fraction) can be reproducibly controlled as 
it is governed by the input parameters. 
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A set of basic operations using specific channel geome- 
tries necessary for manipulating series of bubbles, termed 
"discrete microfiuidics" , has been demonstrated at the 
millimeter scale j5| in a "dry" foam, where the liquid 
content is low compared to the gas content. To down- 
scale these operations and adapt them for microfiuidics 
requires producing a microfoam and information about 
its fiow characteristics. 



Here, we investigate the continuous production of a 
two phase gas-liquid fiow in a fiow-focusing device, and 
the transitions between different regimes of bubble for- 
mation, so as to reach microfoams, thereby extending the 
studies of Ref. |H] to low liquid fractions. In a microfluidic 
fiow-focusing device, a flowing gas thread is forced, by the 
co-flowing surrounding liquid, into a small orifice, where 
the gas thread breaks up at regular time intervals 

At low liquid to total fiow rate ratio this will create a 
microfoam. 



We would like to determine how foam properties (for a 
review, see 0) extrapolate to the micrometer range. We 
investigate here microfoam formation and flow dissipation 
within a microfluidic set-up: a bubble formation orifice, 
followed by a long channel ending with a free exit. We 
finally open perspectives for the study of structure and 
dissipation of ultra-fiat microfoams. 
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2 Materials and methods 

We use a flow focusing geometry: an inlet channel for the 
liquid, another one for the gas, both ending in an ori- 
fice ending up in a straight channel [Tj|3) , 4) t 6), . We produce 
the microfluidic device by soft hthography techniques. We 
first create a mold in a negative photosensitive material 
(SU-8 2100, MicroChem) and then make imprints in a 
polymer (polydimethylsiloxane, PDMS, Sylgard 184) to 
create the actual channel. The PDMS imprint is glued to 
a glass cover slide using a home-built ozone cleaner. The 
exit channel has a height h = 250 fim, vifidth w — 700 
/zm and an orifice width Wor- = 100 ^im. We also use a 
channel where, at a distance of 6 mm after the orifice, 
the foam flows through a second constriction diminishing 
the channel width to 125 fim. We produce the ultra-flat 
foam in a device with a height of about h — 8 /im, width 
w = 400 ^m and an orifice width Wor = 75 /zm; we create 
the mold with a positive photo-resist (ma-P 100, Micro 
Resist Technology). 

For the continuous phase we use deionized water with 
10 % commercial dishwashing detergent (Dreft, Procter & 
Gamble). This solution has a surface tension cr = 38 ± 1 
mN/m, as measured by the Wilhelmy balance method. 
The use of this surfactant resulted in an increased wet- 
tabihty of the solution to the PDMS surfaces [S]. Two 
different syringe pumps were used to push the liquid (11 
Pico Plus, Harvard Apparatus, and KDS 100, KD Scien- 
tific) at flow rates Qi ranging from 4 to 167 /^1/min, with 
±0.5% accuracy. The dispersed gas phase is nitrogen. It 
is driven at constant overpressure Pg (relative to one at- 
mosphere), ranging from 1 to 21 kPa, using a pressur- 
ized tank and a pressure-reduction valve (stability ±0.15 
kPa). Pressure Pin is measured at the entrance of the 
device with piezo-resistive gauge (40PC Honeywell, ±0.2 
kPa accuracy). Since the exit is at atmospheric pressure 
Pout — 1 atm, the overpressure Pg is total gas pressure 
drop Pin — Pout over the microfluidic system (orifice and 
channel) . 

In a typical experiment we vary the gas pressure while 
keeping the liquid flow rate constant. In this way we scan 
the complete pressure range for which bubbles are formed 
at that liquid flow rate. Still images or movies of the result- 
ing flow are then captured with a camera (Marlin F131B, 
Allied Vision Technologies) connected to an inverted op- 
tical microscope (1X70, Olympus), see Fig. ^ for some 
examples. By analysing these we can extract the follow- 
ing bubble quantities: formation rate (break-up frequency 
/), with a precision of a few percent, the bubble volume 
Vb = Abh by measuring the apparent area Ai, occupied 
by the gas in the images, the gas flow rate estimated as 
Qg = Vbf, and the space and time averaged gas velocity 
(u^) estimated as the distance between two consecutive 
bubbles multiplied by /. The edge of the measured area 
Af, (measured with a precision ±1%) is taken in the mid- 
dle of the curved meniscus around the bubble, appearing 
black on images, systematically underestimating (with up 
to 10% for non-touching bubbles) the actual gas volume. 

An important parameter in characterising a foam is 
the foam volume liquid fraction i.e. the proportion of 
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(i) dripping 



(ii) bidisperse 
bubbles 



(iii) bubbly 
fiow 



(iv) alternate 
foam 



(v) bamboo 
foam 



Fig. 1. From wet to dry microfoams. (a) Bubble volume Vj, 
versus gas pressure Pg. Liquid flow rate is kept constant at 
Qi — 167 fil min~^. Numbers identify the different regimes: (i) 
dripping flow; (ii) bidisperse bubbles (two symbols are plot- 
ted for each pressure); (iii) bubbly flow; (iv) alternate foam (2 
rows) with filled symbols; (v) bamboo foam (1 row) with gray 
symbols, (b) Photographs of these regimes. Crossed symbols 
in (a) correspond to pictures in (b). 



the volume occupied by the liquid: 

Vi 



^1 



(1) 



It is estimated by image analysis with ^/ 1 — A^/A. 
Another method is to measure simultaneously the time 
and space averaged bubble velocity {u^) together with the 
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average gas flow rate; the liquid fraction then follows from 

Qg={u^)S[l-^l). 

On the other hand, the proportion of liquid injected in 
the system is 

with Q = Qg + Qi the total two-phase flow. It is sim- 
ply measured as Qi is one of the control parameters. The 
precision in this measurement is at minimum ±10% for 
non-touching bubbles and in general a few percent. The 
two quantities and ai are different because liquid and 
gas can have different velocities. They can be linked using 
the time and space averaged gas and liquid velocities {W^) 
and (ui) since Qi = {ui)S<l>i and Qg = {Ti^)S{l — with 
S the area of the channel cross section. This yields: 



1-ai <Pi 



(3) 



The separate measurement of ai and allows to calcu- 
late the ratio (u7) / (v^) which informs about the relative 
drainage of liquid through the moving foam. The absence 
of relative drainage, (w7)/(wg} = 1, implies that ai = 
while drainage in the direction of the flow, (ul) / (u^) > 1, 
entrains a injected liquid fraction higher than the volume 
liquid fraction a; > 'Pi- 



3 Microfoam formation at low liquid content 



part (dripping flow (i). Fig. ^). At higher Pg, there is 
a coexistence, probably indicating a first-order transition, 
with a second mechanism, where the interface remains in 
the orifice instead of retracting after bubble release. For 
given Pg and Qi, both mechanisms give different volumes 
Vi, V2. This results in a flow of period Ti + T2 0, with 
bidisperse bubbles (ii. Fig. ^). 

Further increasing Pg , we only observe the second pinch- 
off mechanism, always resulting in a monodisperse foam. 
Three possible structures appear, according to the flow 
rate: bubbly flow (iii), alternate foam with two rows (iv), 
or bamboo foam with one row only (v) . No multiple-period 
or chaotic bubbling is observed. This suggests the absence 
of inertial non-linearities during the retraction of the gas- 
liquid interface 'S'. 

For much higher Pg the gas thread stops breaking up 
and a stratified liquid-gas flow is observed. 



3.2 Microbubble volume 

The bubble volume for the foam regimes (iii, iv and v) 
correlates well with the fluid fraction: 



w; 



±_ ^-0.95±0.02 
,3 ' 



(5) 



see Fig. El except for the lowest Qi. At very low liquid con- 
tent ai Qi/Qg <^ \, bubble volumes are approximately 
proportional to [Qi/Qg]^^ ■ 



3.1 Bubbling regimes 

To study bubble formation and the accompanying bub- 
ble topology in the channel we vary the gas pressure at 
constant liquid flow rate. See fig. ^ for examples of the 
observations in the 250 /im high channel, near the orifice 
at the channel entrance for a fixed flow rate of Qi ~ 167 
/il/min. Above a certain threshold in gas pressure Pg bub- 
bles form in the channel. The bubble volume grows when 
increasing Pg, inducing several regimes of bubble forma- 
tion and flow. 

We observe a minimum pressure Pc for which bubbles 
form. For lower pressures the gas-liquid interface does not 
enter the orifice. At this liquid flow rate Pc = 0.9 ± 0.15 
kPa (Fig.^). This effect is probably due to the capillary 
pressure. For a curved interface in the orifice considering 
the limit of bubble detachment, the Laplace pressure of a 
wetting interface is 

Pa = cr ( — + —)= 1.1 ±0.2 kPa, (4) 
\ri 7-2/ 

where ri = h/2 and r2 — Worj'^ are the principal radii of 
curvature, of the same order as Pc- 

Above the initial pressure Pc, a gas thread is forced 
into the orifice and fills a bubble after the orifice. This 
thread pinches off and releases the bubble. After break- 
up, the gas-liquid interface retracts to its initial position, 
as reported in j4j , returning completely into the upstream 




Fig. 2. The bubble volume 14, in units of ni'o^. — 1 nl, depends 
only on the injected liquid fraction ai. The solid line is a linear 
fit to all data except the lowest Qi: log(y/Wo,,) = (—0.95 ± 
0.02) logo, + (1.78 ±0.05). 

This correlation is similar to the one observed in ax- 
isymmetric conditions by (QulOj who measured T4/Wor ^ 
(Oi/Qg)^^'^^^" *'^, in the opposite case of high liquid frac- 
tion with Qi/Qg > 5, with separated bubbles, while here 
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bubbles are in contact or in short spacing after forma- 
tion. It differs from the Vb ^ Pg/Qi scaUng observed by 
Garstecki et aL ^ : we will see below that the gas flow rate 
and applied pressure are not proportional for a low liquid 
fraction microfoam in a channel. Moreover ref. j3j inves- 
tigates high liquid fraction flows at channel aspect ratio 
of 0.04 (flatter channels) much lower than the present 0.3. 
Both points could explain the experimentally observed dif- 
ferences. 

In this geometry, one or two rows of bubbles are ob- 
served. The transition from bubbly flow to a foam in which 
bubbles touch each other is governed by the volume liq- 
uid fraction which is related to ai through equation 13 
Therefore increasing the number of bubble rows when the 
foam state appears would either require decreasing Vb at 
given (Pi or enlarge the space for newly formed bubbles. 
The former can be achieved for instance with an orifice as- 
pect ratio Wor/h closer to 1 to restrict liquid flow; and the 
latter by a lower ratio of orifice to channel width Wor/w. 



can be used to define a non-dimensional frequency, known 
as the Strouhal number. 



St- 



/ 



(7) 



After rescaling by // all data collapse on a single curve 
where St is a function of gas pressure only. Two regimes 
are observed: for low gas pressures in the case of bubbly 
flow we find St ~ 0.1(Pg — Pc)/{2a/wor), while for higher 
Pg, the Strouhal number saturates to a constant value of 
St = 0.16. 

We infer that these two regimes are the consequence 
of two stages during the bubble formation: 

1. Gas filling of the orifice: At low Pg the frequency varies 
like / ~ MPg - Pc) ~ Qi{Pg - Pc). In other words, 
the period is proportional to a characteristic time that 
varies as 



T 



mpg-Pc)]- 



(8) 
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Fig. 3. The bubble formation frequency / in units of /; = 



gas pressure Pg in units of 2a/wor- It shows 



plateau for high Pg-. the dotted line is a piecewise linear fit. 
Symbols correspond to different values for the liquid flow rate 
Qi: (circle) 4, (diamond) 15, (square) 20, (right-pointing trian- 
gle) 30 and (down-pointing triangle) 40 /il min~^, for bubbly 
flows (open symbols), alternate foams (filled symbols), bamboo 
foams (gray symbols). 



3.3 Microbubble formation frequency 

The frequency / of bubble formation (Fig. |2Jl first in- 
creases linearly then reaches a plateau for increasing Pg. 
The typical time and frequency linked to the liquid fiow 



We interpret this time as the time necessary for the gas 
to fill the orifice, prior to break-up. It decreases with 
increasing Pg — Pc since the gas pushes the fluid with a 
velocity increasing with pressure. Tg also decreases for 
increasing Qi with the flow focusing conflning more 
and more the available space for the gas thread. Note 
that the relation / QiPg was proposed by 0, ver- 
ified for varying Qi but with a constant Pg-. here we 
also investigate the effect of gas pressure. 
Liquid mediated thread pinch-off: For high Pg the bub- 
ble formation frequency is proportional to / ~ Qj. The 
period only depends on 



Tl 



Ql 



(9) 



the time to pinch off the gaseous thread when the liq- 
uid flow is blocked by the bubble at the outlet 4,. 

The bubbling period at low pressure is limited by the 
gas flUing, while it is limited by the liquid driven thread 
contraction in the high gas pressure regime. Note that the 
transition from a Tg to a ti dominated break-up frequency 
is accompanied by the regime change from bubbly flow to 
foam. In the foam state, the liquid flow restriction seems 
more efficient (see figure iv and v). 

As a conclusion, there are two stages during formation: 
the first associated with the filling of the orifice by the 
gas (Tg) and the second reflecting the pinch-off of the gas 
thread (t;). They have different gas pressure dependency 
(Tg depends on gas pressure while r; does not), which cre- 
ates a cross-over apparent in the bubbling period 1// that 
depends on Tg and t;. 



4 Foam flow 



(6) 



We now turn to the flow of a foam in the microchannel 
after formation. When we measure the average gas flow 
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rate as a function of the apphed gas pressure drop Pg, we 
observe a highly non-hnear response (see Fig.0J. We find: 
a threshold (section |4.1() . a non-linear slope (section |4.2|I 
and a discontinuity upon the transition from an alternate 
to a bamboo foam fsection f4.3|l . The pulsation of the flow 
rate at the bubble formation frequency is finally examined 
(section I4.4|l . 



4.1 Pressure threshold 

Fig. 2] shows a threshold in pressure for the establishment 
of bubbly flow. It is found to be 1.0 ± 0.1 kPa, if this 
parameter is left free in the fit. It is compatible with the 
above explanation by a capillary effect (1.1 ± 0.2 kPa 
according to equation 0} at the orifice. 

In presence of a second constriction (data not shown), 
we obtain for regimes (iii) and (iv) the same result as in 
Fig. Q), translated by about 0.45 kPa along the P-axis 
(compatible with the expected Laplace pressure necessary 
to overcome the second constriction, 0.6 kPa). This con- 
firms that the threshold is induced by capillary effects. 
On the other hand, the slope originates from dissipative 
effects in the channel. We thus write the total pressure 
drop as the sum of two contributions, 



channel i 



(10) 



where is the static orifice contribution and APchannei 
is due to dynamic dissipation in the channel. 



300 




2 3 4 5 6 

gas pressure [kPa] 

Fig. 4. Gas flow rate Qg vs. gas pressure Pg showing a cap- 
illary threshold and discontinuity at the transition from alter- 
nate to bamboo foam {Qi = 20 fA min~^). 
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Fig. 5. Graph using adimensional axes Ca = ^Qg/Sa ver- 
sus APchannei = Pg — Pc- Data for different liquid flow rates 
are superimposed, Qi = (o) 4, (□) 20, (>) 30 and (v) 40 /A 
min~^. The solid lines on both plots are flt to data on alter- 
nate and bamboo foams at Qi — 20 fi\ min~^ with the power 
law APchannei /(cr/h) = PCa^^^, with f3 = 1.7x 10^ (alternate) 
and /3 = 2.5 X 10^ (bamboo). 



4.2 Non-linear flow-rate to pressure dependence 

The flow rate is highly non-linear above the threshold, 
see figE] The gas flow rate increases faster than a linear 
function of pressure for the alternate foam structure. 

It can be interpreted by assuming that dissipation mainly 
occurs in the liquid films, close to the walls [2]. Introduc- 
ing the capillary number Ca = fiv/a (of order 10~^— 10^"^) 
containing the bubble velocity v (estimated as w ~ Qg/ S 
in the dry foam state) and liquid viscosity /i, the pressure 
drop writes 



^-Tchannel — A — aUa ' , 



(11) 



with n the total number of bubbles, Lproj the projection 
on the cross section of the wetting perimeter per bubble, 
S the cross section area, and A a numerical constant jll|. 
The effect of the orientation of the films relative to the 
foam movement is included in the Lproj variable {Lproj = 
L cos a for a film whose normal vector is slanted by an 
angle a with respect to the flow direction). 

In each regime, a fit of the pressure drop APchannei by 
Ca}l^ gives a correct agreement, see figure Compared 
to a Newtonian flow whose drag pressure grows propor- 
tionally to the flow rate, the drag pressure grows with a 
lower exponent of the flow rate. It is an effect of the lubri- 
cation films between bubbles and walls that thicken when 
flow rates increase, a phenomenon described first in |12| . 
for the motion of a single in a capillary tube. 

We measure the projected friction length Lproj as the 
projection of the length between vertices centers, for con- 
tacts with the walls parallel to the image plane, neglecting 
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side walls. We deduce for the numerical constant the value 
A = 22 ± 5. The liquid coflow, varying over a decade, does 
not influence much the gas pressure drop in our geometry 
and the value of this constant (see figEJ. 

We can compare the constant A between our microchan- 
nel and the millimetric channels with a comparable as- 
pect ratio (but without liquid coflow) studied by , who 
found A = 38 ± 04. We observe less friction: in the present 
experiments liquid is injected continuously, possibly ex- 
plaining wetter foams and a lower value for A. Since we 
see no visible change in the thickness of lubricating liq- 
uid films between gas and walls, which would change drag 
forces, we assume that liquid flows mostly in the corners 
\n\. Thanks to this corners the liquid flow is not as ob- 
structed by bubbles in the channel as with a cylindrical 
geometry [T^ . 

We find that the measured foam liquid fractions is close 
to the injected fraction <l>i ~ a;: the liquid fraction is here 
controlled by the injection parameters. Using (O we ob- 
tain (ui) I {u^) ~ 1: there is no substantial relative drainage 
of the liquid within these foams. From conservation of liq- 
uid mass in the corners, we deduce that an increase of 
the liquid flow rate with no variation in average velocity 
should lead to an increase of the diameter of the corner. 



micrometer scale. The forcing frequency is therefore grad- 
ually attenuated, allowing to consider the flow as steady 
further in the channel, as if pushed by a steady pressure. 
In view of the use of foams in channels with obstacles, 
this would prevent resonance effects with the frequency of 
passage of bubbles on the obstacles downstream. 




Fig. 6. Space-time diagrams of the foam flow in the constric- 
tion sample: (a) at the channel entrance; and (b) just before the 
exit constriction. The vertical axis is the time, flowing down- 
wards; the horizontal axis is a (small zone of) the axis of the 
channel, with the foam flowing from left to right, Qi = 16.7 \A 
min~^ and P — 4.9 kPa. Dark pixels indicate a bubble edge. 



4.3 Discontinuities in the flow-rate 

The structure transition from an alternate to a bamboo 
foam induces a discontinuous decrease of the gas flow rate 
(Fig. II. 

This is a signature of the discrete character of the 
foam: We find by image analysis that the transition to 
bamboo structure is associated with an increase in the 
bubble number n of 15 % (compactification). It is also 
associated with an increase of the projected length Lproj 
of 50 %. Both factors are consistent with the observed 
increase of APchannei by 50 %, considering equation ^2 
that states that APchannei is proportional to nLproj- We 
conclude that the rearrangement of bubbles induces the 
discontinuity in the pressure drop. 

A consequence of this findings is that, in some ranges, 
a given flow rate (here for Qg in between 170 and 260 
/Lil min~^) can lead to two possible foam states, each one 
associated with a different pressure drop. 

4.4 Flow pulsation at bubbling frequency 

The above flow rate measurements are in fact time-averaged 
flows. The foam velocity at entrance oscillates at frequency 
/ (Fig-Et)j between 3.1 and 8.3 cm/s. However, this oscil- 
lation is damped over a few mm along the channel (data 
not shown) . Very little oscillation is observed at the chan- 
nel exit (Fig. Eb). 

Due to the periodic creation of bubbles, the pressure in 
the entrance orifice probably varies at frequency /. Each 
new bubble has to push the foam to create a place for it- 
self; thus inducing the flow oscillation. Friction, probably 
mainly at the channels walls, homogenizes the flow at the 



5 Ultraflat foam: distortion effects 

In order to downscale even more flowing foams, we reduced 
the channel height to produce ultra- flat foams. The ultra- 
flat channel presents a 30-fold decrease in height, to 8 
/im, and a 18-fold decrease in aspect ratio for the channel 
section, to 0.02. We continuously produce various foams 
in such a channel (Fig. 0, including one with 3 bubble 
rows (Fig. [7^). 

We can dry it, using the following batch method. We 
shut the liquid inlet and pull the syringe at the gas inlet. 
As long as the underpressure is smaller than the Laplace 
threshold (here 11 kPa), the bubbles are blocked by the 
orifice and only liquid flows out of the foam. This forced 
drainage yields hexagons with a standard deviation in the 
edge length of only 1.8 % (Fig. Since the apparent 
wall thickness on images (10 /im) is comparable to the 
height, the bubble walls are probably very curved (when 
looking at their profile on a cross-section perpendicular to 
the image plane), and with no flat film between bubbles, 
contrary to the previous set-up. Thus the actual liquid 
fraction of the central bubble row is probably smaller 
than, but close to, the apparent one (fraction of black 
pixels) 'P « (pY^°-9& ^ 10~i. This contrasts with usual 
foams with larger aspect ratio, where the same picture 
of hexagons with straight walls and small vertices would 
correspond to much lower fluid fractions, <Pi ~ 10^^ < 

Foams flowing in this ultraflat channel undergo an un- 
usual boomerang-like distortion, with films near the side 
walls pointing backwards (Fig.|3:,d), the opposite of foams 
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with no coflowing hquid |S] • We indeed expect a larger fric- 
tion at the center than on the edges, where bubble walls 
along the top and bottom plates are thin, than on the 
channel sides, where water accumulates. These edges are 
likely more entrained by a faster liquid flow than the cen- 
ters {{ui)/{u^) > 1, assuming that for small aspect ratios 
corner sections are less expandable). Within the reference 
frame moving at the liquid velocity, the bubbles move in 
the opposite direction, from right to left. The distortion of 
bubble shapes in this reference frame, with films near the 
walls pointing forward, is then similar to the observations 
of millimetric bubbles pushed without any liquid flow 5 . 




Fig. 7. Ultraflat foams in the 8 ^m high channel: (a) flowing, 
wet and (b) static, dry 3-rows foams; (c) 2-rows and (d) 1-row 
flowing dry boomerang foams. Flows from left to right. 



6 Conclusions 

We describe the formation and flow of a foam in a conflned 
microchannel. The transition from bubbly flow to foam de- 
pends on $1 which is governed by the interplay between 
control parameters Pg and Qi and the channel geometry. 
The frequency at which bubbles are formed behaves dif- 
ferently for bubbly flows and for foams. For foams the 
formation frequency only depends on the liquid flow rate 
controlling the speed at which the gas thread is pinched 
off. Foam flow and bubbly flow in microchannels are highly 
non-linear. The flow focusing oriflce induces a threshold Pc 
due to capillary effects in the flow-rate to pressure char- 
acteristic. The data for both bubbly and foam flow give 
good agreement to Pg ~ Pc Ca?'^^. The prefactor in this 
relation depends on the dissipation in the channel related 
to the topology. 

Microfoams can be down-scaled to as small heights as 
8 fim.. The liquid fraction can be varied continuously over 
the complete range from the dry to the wet limit. We see 
an unusual deformation for foams flowing in this channel 
probably caused by relative drainage. 

The interplay between geometrical parameters (chan- 
nel aspect ratio, oriflce aspect ratio and oriflce to channel 



ratio) merits more attention. It governs various effects as 
the transition from bubbly flow to foam, the foam topol- 
ogy, and thereby dissipation in the channel, and the dis- 
tortion of the foam cells. Comprehending this interplay 
will be a necessary step in the development of microflu- 
idic applications in which the straight channel section will 
be replaced by more complex geometries allowing opera- 
tions like mixing, separation, breaking and coalescence of 
bubbles. 
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